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Quasiclassical theory of disordered Rashba superconductors 
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We derive the quasiclassical equations that describe two-dimensional superconductors with a 
large Rashba spin-orbit coupling and in the presence of impurities. These equations account for the 
helical phase induced by an in-plane magnetic field, with a superconducting order parameter that 
is spatially modulated along a direction perpendicular to the field. We also derive the generalized 
Ginzburg-Landau functional, which includes a linear-in-gradient term corresponding to the helical 
phase. This theory paves the way for studies of the proximity effect in two-dimensional electron 
gases with large spin-orbit coupling. 
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I. INTRODUCTION 

Breaking time-reversal symmetry in superconductors 
that lack inversion symmetry can result in new and 
interesting phenomena, such as topological phased or 
magneto-electric couplingj^. For instance, the ap¬ 
plication of an in-plane magnetic field to a non- 
centrosymmetric superconducting layer allows for the for¬ 
mation of a “helical” phase, which is characterized by 
a spontaneously modulated order parameter in the di¬ 
rection transverse to the fielcP. On the phenomenolog¬ 
ical level, this helical phase results from the fact that 
the Ginzburg-Landau functional may admit a linear-in- 
gradient terirP. This term, also known as a Lifshitz in¬ 
variant, describes the coupling between the magnetic field 
and the supercurrent in the presence of spin-orbit cou¬ 
pling. Thus, the helical state is conceptually different 
from the no nun iform Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) statd^, which arises from a change of sign of 
the quadratic-in-gradient terirPl that fixes the amplitude 
of the modulation, but not its direction. 

The amplitude of the Lifshitz invariant was calculated 
for the microscopic model of a two-dimensional simercon- 
ductor with Rashba coupling in the clean limiiP. This 
study showed that the modulation wavevector in the he¬ 
lical phase is proportional to the spin-orbit-induced dif¬ 
ference in the densities of state of electrons with opposite 
helicities. The parallel upper critical field in the helical 
phase was computed both in the cleanand diffusiv^^i^ 
regimes. The relevan ce of this effect for superconductiv¬ 
ity at oxide interfaceJ^^Ull and in single-atomic layers of 
PtP 

was also discussed. Furthermore, the possibility of 
a helical phase in three-dimensional superconductors was 

debatecP^ltn!_ 

Here we provide a microscopic theory of the helical 
phase in the disordered case. While homogeneous super¬ 
conductivity is basically unaffected by disorder as long as 
the localization length remains m uch la rger than the co¬ 
herence length (Anderson theorerrfiSlIIl) ^ the FFLO state 
is rapidly destroyed by even weak disorder!^. By con¬ 
trast, the helical phase is modified by the disorder, but it 
survives up to fairly large disorder strengths. The qua¬ 


siclassical theory offers a convenient framework to dis¬ 
cuss the equilibrium and out-of-equilibrium properties 
of superconductor^. The difficulty to elaborate such a 
theory in the case of non-centrosymmetric superconduc¬ 
tors is related with the large spin-orbit induced splitting, 
which is responsible for the helical phase. Namely, the 
normal-state dispersion of the material consists of two 
so-called helical bands with different densities of state at 
the Fermi level. In the clean limit, a large spin-orbit cou¬ 
pling results in an effective two-band superconductivity 
with Cooper pairs that form separately in each of the he¬ 
lical bands, while sharing the same gap functiorP^. The 
paramagnetic effect induced by the magnetic field results 
in a modul ation of the order parameter similar to the 
FFLO phas(P^. The preferred modulation wavevec¬ 
tor is opposite in the two bands. Thus, the modulation 
results from a competition between the two bands and 
their difference in densities of states is crucial in order to 
obtain a non-vanishing result. Quasiclassical Eilenberger 
equations have been derived to describe this situatiorPSl. 
Their solution permitted to consider the crossover be¬ 
tween FFLO and helical phases as well as to predict the 
density of states in the helical phase. However, impurity 
scattering mixes the bands. Thus, inter-band super¬ 
conducting correlations may develop. So far quasiclas¬ 
sical equations for the diffusive regime were deri ved in 
the spin basis but neglecting the spin-orbit splittin^^iissj 
Thus, they did not account for the helical phase. More 
recent works considered the case of very weak spin-orbit 
couplin^2ll27l^ where only a higher order effect is expected 
to survive. 

In our work, we derive quasiclassical equations in the 
helical basis which account for the spin-orbit splitting 
and therefore describe the helical phase. To set the stage, 
we first derive the quasiclassical equations in the helical 
basis when the spin-orbit coupling is small and the helical 
modulation may be neglected. Within that approxima¬ 
tion, we demonstrate that for moderate disorder, despite 
inter-band impurity scattering. Cooper pairing predomi¬ 
nantly occurs in separate bands. Thus, we may derive a 
set of quasiclassical equations for each helical band sepa¬ 
rately. This allows us to take into account the spin-orbit 
splitting and obtain the helical phase. The equations 
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we derive describe the entire superconducting phase at 
arbitrary temperatures and magnetic fields. We also de¬ 
rive the dependence on disorder of the Lifshitz invariant 
in the generalized Ginzburg-Landau functional to further 
characterize the transition to the helical phase. 

II. MODEL 

A two-dimensional electron gas with Rashba spin-orbit 
coupling is described by the Hamiltonian 

7/o = ^ Kk + a(k X z) • cr] ak. (1) 

k 

Here, Ok = are annihilation operators for 

an electron with spin a =t,4- and momentum k, ^k = 
k^/(2m) — fi is the kinetic energy associated with the ef¬ 
fective mass m and chemical potential fJ., a = Ago/ y/2mfi 
characterizes the spin-orbit coupling strength, and z is 
a unit vector normal to the electron gas. Furthermore, 
Pauli matrices ai act in spin space. 

The Hamiltonian Q can be diagonalized in the so- 
called helical basis, 

77o = ^ 4 (?k + aka^) Ck, (2) 

k 

where Ck = U-^au, with Uu = ^ ^ . 

Here, (fk is the angle of k in the (a;?/)-plane. According to 
Eq. (i), the helical bands are characterized by the same 
Fermi velocity, v = y/2^/m + and different densities 
of states, v\ = u(l — Xa/v), where u = ml (27r) and A = ± 
labels the helical bands. 

A magnetic field B acting on the electron spin induces 
a Zeeman term, 

Hz = aj^h • crak = ^ cj^hk ■ crck, (3) 

k k 

in the spin and helical basis, respectively. Here, h = 
and hk = (h^, —hy • k, —(hy x z) • k)^, where hy 
and hz are the in-plane and out-of-plane components of 
h, respectively, while k = k/A:. 

Superconductivity is described by a BCS Hamiltonian 
with s-wave pairing^^, 

Hs = -Es «kt«U+qt«-k'+qfak't, (4) 

k.k',q 

where Vs is the pairing constant. Within mean field the- 
ory, Eq. @ can be expressed in the helical basis, 

Hs = \Y^ ^A(’^)^(q)4Aclk+qA + h-c., (5) 

k,q,A 

where tA(k) = while the order parameter solves 

the self-consistency equation 

^(q) = - ^Vsy^,tx{k)* (c-k+qxckx) ■ (6) 

kA 


Finally, potential scattering by impurities is accounted 
for by the term 

i?imp = ^a|,Ho(k-k')ak' = ^c|,C/o(k-k')W'k.k'Ck', 

k,k' k,k^ 

(7) 

where Hk.k' = UlUk' and the disorder potential satisfies 
the average (t/o(q)t/o(-q'))dis. = C7^<5q,q'- 

Introducing annihilation operators CkA = 
(ckAj AA(k)c^i^^)'^ in particle-hole space, the total 
Hamiltonian takes the form H = ^ C^Hk.k'Ck' 

with Ck = (Ck+Xk-)^ the Bogoliubov-de Gennes 
Hamiltonian 

■Hk.k' = [('^k + ak<Tz) Tz + h.k- cr] (5k,k' (8) 

+t7o(k - k')IFk.k'T. + A(k - k'). 

Here, A(k - k') = A(k - k')T+ -h A*(k - k')T_ and r, 
are Pauli matrices acting in particle-hole space with t± = 
(ja: ± iTy)/2. 

We now define the Matsubara Green’s function Q 
which solves the equation {ioj—'H)G = 1, where w = (2n-|- 
l)7rr (n integer) is a Matsubara frequency at tempera¬ 
ture T. Treating the disorder within the self-consistent 
Born approximation allows one to relate the impurity- 
induced self-energy to the disorder-averaged Green func¬ 
tion G = {G)dis.- In a mixed momentum/space represen¬ 
tation, this relation reads 

E™P(r) = C/oV,^Wk.k'Gk'(r)WkykT., (9) 

k' 

provided that the Fermi wavelength is much smaller than 
the scale for spatial variations of Gk(r). The gap equa¬ 
tion ([^ yields 

A(r) = -^PsT5]Tr[r_Gk(r)]. (10) 

cj,k 

The Dyson equation for the Green function Gk(r) can 
be simplified within a quasiclassical approximation that 
assumes a separation between the large energy scales 
characterizing the electron spectrum in the normal state 
and the small energy scales associated with the supercon¬ 
ducting properties. We will first derive the quasiclassi¬ 
cal equations for the full matrix Green function. Impor¬ 
tantly this requires to neglect the difference between the 
densities of states in the two helical bands. Thus, these 
equations are valid only in the case of a small spin-orbit 
coupling, a V, but at arbitrary disorder. We will then 
discuss how to take into account the difference between 
the densities of states to describe the helical phase at 
strong spin-orbit coupling, a < v, and moderate disor¬ 
der, 1/r < Aso- 

III. WEAK SPIN-ORBIT COUPLING 

In the regime of weak spin-orbit coupling, a <C u, the 
splitting of the Fermi surface between different helical 
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bands is negligible. In that case, the derivation of the 
quasiclassical Eilenberger equations solved by the energy- 
averaged Green function, 


5 (r,k) = ^ / (i4 T^Gk(r), 


( 11 ) 


is standard. We obtain 


- uk • Vg(r, k) = [(w -I- iA(r) -|- ihk ■ cr) + iAgoCTz 


+a“P(r,k),g(r,k) 


( 12 ) 


with the reduced self-energy 


a™P(r,k) = ^C/l(C7k-g(r,kOC/l,)fc,C/k. (13) 


Here 1/r = 27^11^ and (...)£ denotes averaging over the 
Fermi surface. Using the explicit form of the matrices C/k, 
one may write Eq. (131 as 


cr™P(r,k) = — |(g(r,kO)k/+c^:r(5(i',k0)k'O'^ (14) 


-k- 


-f izx)[(k'5r(r,kO)k/iCrx]) }• 


The self-consistency equation (10) reads 
ittYsv r 


A(r) = 


T^Tr r_(g(r,k))i. 


(15) 


Moreover, the normalization condition g^(r, k) = 1 must 
be satisfied. 

At large disorder, l/r ^ h,A,T, the quasiclassical 
Green function is essentially isotropic over the Fermi sur¬ 
face. That is, we may approximate 5 (r,k) « g(r) -|- k • 
g(r), where the second term is a small anisotropic correc¬ 
tion. The normalization condition imposes = 1 and 
{g,g} = 0. If, furthermore, we assume Ago ^ h,A,T, 
Eq. (12) averaged over the Fermi surface yields 


iAso(Jz + ^{9 + crxg(Jx),g 
4t 


= 0, 


(16) 


Using {g,g} = 0, we find that Eq. 0 is solved with 

1 


g = -vTgoVgo + 2 zt ( 1 + 

hn ^ 1 


-I- 


2rAs 


A^2A2 ; (h|| X z)(T,5o[A,5o] 

^ ‘ ^SO 

(^v9o] [Tz,go]- (18) 


As a second step, we find the isotropic correction Sg 
by considering the terms in the next-order correction to 
Eq. (16) that are non-diagonal in helical space. We ob¬ 
tain 


6g = -i 


2Ag 


1 


2tAb 


<Jx9o] [Tz,go]- (19) 


Finally, taking into account Eqs. (18) and (19) to obtain 


the diagonal terms, we find that gQ obeys a standard 
Usadel equation, 


DV ■ (go^go) = 


(w -f iA)Tz + -TzgoTz, go 


where D = iir/m is the diffusion constant and 


F = 2Thi 1 -h 


2t2A2 


-f 




2rA2 


( 20 ) 


( 21 ) 


is an effective pair-breaking parameter. Note that the or¬ 
bital effect of the magnetic field may be readily accounted 


for in Eq. (20) with the substitution V —>■ V — *eA[T 2 ,.], 


where A is the vector potential. 

Properties of superconductors described by Eq. (20) 
are well establishecP^. In the absence of a magnetic field 
F = 0, and therefore the superconducting properties are 
insensitive to the spin-orbit coupling. A finite magnetic 
field suppresses superconductivity. At T = 0, a second- 
order transition to the normal state occurs at the upper 
critical field determined by the equation F-l-Forb = Ao/2, 
where Forb = eDBz and Aq is the superconducting gap 
at zero temperature in the absence of the magnetic field. 
Thus, the in-plane critical field is given as 


h\\r 


Ag, 



in leading order in 1/t and Ago. Thus, it is satisfied if 
g{r) = goir) + Sg{r), with [a^,go{r)] = [cr 2 , 5 o(r)] = 0 
and dg(r) < goir). 

We therefore look for a solution go{r) that is diagonal 
in helical space and obeys the normalization conditions 
gl = I and { 50 , 9 } = { 50 , < 55 } = 0- As a first step, 
we determine g by inserting the expansion for g(r,k) 
into Eq. (12) and averaging over the Fermi surface after 


multiplication with k. This yields 


-vVgo - [{ax -I-fzx)h||cr2T2,5o] 


(17) 


= [iAsodz + ^go,g] - ^[idx + iiy-)[g,dx\,go\- 


In the out-of-plane configuration, the paramagnetic 
mechanism yields a critical field = Aso\/A qt, which 
is much larger than the in-plane critical field /ly,, for 
1/r ^ Ago. However, typically the out-of-plane criti¬ 
cal field is determined by the orbital mechanism yielding 
^±c ~ |5/^B</o/(27r^^), where (po is the flux quantum and 
is the superconducting coherence length. 

Using these results, we can now check our assump¬ 
tions used in the derivation of the Usadel equation. We 
find that the approximations g, Sg ^ go are well obeyed 
provided that Aq ^ 1/r A^o/Aq. To describe the 
system at larger disorder, 1/r A^^/A, the spin ba- 

si^^is better suited than the helical basis. In the regime 
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Ago ^ l/r ^ Ag^/Ao, both approaches are valid and 
yield the same results. 


Further examining Eq. (18l, we note that, at 1/t ^ 


Ago, the anisotropic correction to g is essentially diagonal 
in helical space. This suggests that, in this regime, we 
may project the Green function onto the helical bands, 
thereby neglecting the small inter-band superconducting 
correlations. In doing so, we then are able to account for 
both the large intra-band superconducting correlations 
and the difference in the densities of states between the 
helical bands, and thus access the helical phase at strong 
spin-orbit coupling, a ^ v. 


IV. STRONG SPIN-ORBIT COUPLING 

As shown above, at “moderate” disorder, 1/r <C Ago, 
superconductivity takes place separately in each of the 
helical bands. In that case, we may treat the spin-orbit 
splitting non-perturbatively by deriving the quasiclassi- 
cal equations for the Green functions projected onto the 
helical bands. 


5 ( 1 ', k) 


5 +(r,k) 0 

0 5-(i’,k) 


(23) 


where gv (r, k) are matrices in particle-hole space only. 
Inserting the ansatz (231 into Eq. and performing the 
quasiclassical approximation, we find 


-vk-dxgA(r,k)= (uj + iA(r)) + a‘^^(r,k), gx(r,k) 

(24) 

where c?a = V — i(A/v)(hji x z)[tz, .] and 


A' 


(25) 


+AA'k-(k'gA'(r,kO)k' • 


Here, the different densities of states in the helical 
bands result in different scattering rates 1/tx = 2niy\UQ. 
The different densities of states also appear in the self- 
consistency equation (101 that now reads 


A(r) = 


iirVs, 


r^j^Atr [r_(5A(r,k))k 


(26) 


As usual, the normalization condition 5 ^(r,k) = 1 must 
be satisfied. 

In the presence of disorder, we proceed as in the case 
of weak spin-orbit coupling. Namely, we approximate 
5A(r,k) « 9\{r) + k ■ gxir) and gxir) « go{r) + X6g{r), 


with go = 1 and {go,Sg} = {go,g\} = 0. Inserting this 
expansion into Eq. ( |24[ ) and performing the averages over 
k as discussed above, we find 

g± = ( — 9±5o ±* — (^11 X z)[Tz,go] ), (27) 

yr± UTzp / 


Sg = -—goid+-g+-d--g-), 


(28) 


where 1/t l/(2rA). Finally, we obtain the Us- 

adel equation 


Dd ■ {godgo) = 


{u} -h iA)rz -k -T^goTz, go 


(29) 


where D = t{v'^ + c?)l2, 5 = V -I- 2fa/ (c? -I- r'^)(/i|| x 
z)[tz, .], and f = 2r/iy (u^ — a^)/(r!^ -I- a^). Equation ([^ 
which describes Rashba superconductors in the regime 
A ^ 1/t ^ Ago is the main result of our work. 

In addition to a pair-breaking parameter, F, which re¬ 
duces to Eq. (21) in the limit a v and 1/t ^ Ago, 
the Usadel equation (29) contains an effective vector po¬ 


tential. This vector potential may be eliminated by the 
unitary transformation go(r) = U(r)go(r)H^(r), where 
U(r) = exp[— (f/2)q • r Tz] and 


q= - 


4a 


r(^|| X z). 


(30) 


Under this transformation, the order parameter trans¬ 
forms as A(r) —>■ A(r)exp[— fq- r] = const, which cor¬ 
responds to a spatial modulation. This is the so-called 
helical phase. The modulation is necessary to realize the 
ground state with zero currenllH]. Furthermore, it leads 
to an increase in the upper critical field 


^(q) / Ao + 0,2 


(31) 


Note that the out-of plane component of the magnetic 
field does not contribute after projecting the Green func¬ 
tions onto the helical band^^ . 


V. GINZBURG-LANDAU FUNCTIONAL 

Near the second-order transition to the normal state, 
the order parameter vanishes. In that regime, it is suf¬ 
ficient to solve the Eilenberger equation (24) perturba- 
tively in A around the normal state solution, g^^^ (r, k) = 


TzSign(w). Inserting the perturbative solution up to the 
third order into the gap equation (26), we obtain as the 


result of a long, but straightforward calculation 
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1 


A(q) = 2t:T 


E 

a ;>0 


1 v'^q^ / 1 \ 

uj uj^{uj + ri) 8w^f2 \ 2t{uj + fi)y 


a 

u}^{uj + n) 


q- (h|| X z) 


A(q) 


1 


E A(qi)A(q2)A*(q3) 


qi+q2-q3=q 


(32) 


assuming h,vq <C IIt,T <C Ago- Here we used notation H = w + l/(2r). For temperatures close to the critical 
temperature, this equation minimizes the free energy functional 

T = i.Jid^r) I + 7/if) |Ap + |A|" + /3|VA|2 - (hy x z) (A*VA - AVA*)| , (33) 


q (t) / qo 



FIG. 1: Helical modulation wavevector [in units of qo = 
2q:/i||/u^] as a function of the disorder strength, for various 
amplitudes of the spin-orbit coupling, near the superconduct¬ 
ing critical temperature. 


where T^o is the critical temperature at zero magnetic 
field, P = (^TuV4) E.>o[1 + + 4a;r))]/(u;2H), 

and K = aj = 2a7rT^^^Q l/[a;^ (w-I-H)]. The mod¬ 
ulation of the helical phase is easily obtained as q = 
— H;/(2/?)h|| X z; its dependence on the disorder strength 
is shown in Fig. In particular, this result confirms 
that, in the lowest order in a/v, the helical modulat ion 
wavevector only varies by a factor 2 between the clearP^^ 
and dirty limitflil. Note that the critical temperature is 
shifted downward, Tc{h\\) = Tco(l — y/iy -h Pq'^). At any 
disorder, the shift vanishes when a —>■ u (a result that 
was found previously in the clean cas^^^ , while the heli¬ 
cal modulation wavevector remains constant. 


VI. CONCLUSIONS 

In this work, we showed that the helical phase that 
arises in two-dimensional Rashba superconductors can 
be described with quasiclassical equations, both in the 
clean and moderately disordered case. Importantly, the 
difference between the densities of states in the helical 
bands can be accounted for within quasiclassics because, 
even though impurities lead to inter-band scattering, 
superconductivity remains essentially intra-band when 
the disorder-induced broadening of the helical bands does 
not exceed their energy splitting. 

Beyond the bulk properties discussed here, a natural 
application of this formalism will be to investigate the 
proximity effect between a conventional superconductor 
and two-dimensional semiconductors with a large spin- 
orbit coupling such as, for example, InAs. Furthermore, 
it will be of interest to generalize the formalism to differ¬ 
ent types of spin-orbit coupling as well as to the three- 
dimensional case. 
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